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Abstract 

This work is concerned with the quasineutral limit of the Vlasov-Poisson system in two 
and three dimensions. We justify the formal limit for very small but rough perturbations of 
analytic initial data, generalizing the results of [12) to higher dimension. 


1 Introduction 


In a non relativistic setting the dynamics of electrons in a plasma with heavy ions uniformly 
distributed in space is described by the Vlasov-Poisson system. Throughout this paper, we will 
focus on the 2 and 3 dimensional periodic (in space) case. We introduce the distribution function 
of the electrons for t € M + , ( x,v ) € T d x M. d where is the d-dimensional torus and 

d = 2,3. As usual, f(t,x,v)dxdv can be interpreted as the probability of finding particles with 
position and velocity close to the point (x, v ) in the phase space at time t. We also define the 
electric potential U(t,x ) and the associated electric field E(t,x). 

We introduce the positive parameter e defined as the ratio of the Debye length of the plasma 
to the size of the domain. The Debye length can be interpreted as the typical length below 
which charge separation occurs; it plays an important role in plasma physics. For a more 
detailed discussion on this subject we refer to the introduction of [ID] . Adding a subscript in 
order to emphasize on the dependance on e. we end up with the rescaled Vlasov-Poisson system: 


dtfe + V • V x f E + E e • V v f e = 0, 

E e = — V x U e , 

£ = f^ d f e dv / T d X Rd /e dv dx, 

fe\t=0 = fo,e — 0, ffd x ^dfo,edxdv = 1, 


( 1 . 1 ) 


and the energy of this system is 

£{fe{t)) ■= \ [ f e \v\ 2 dvdx + %- f \V x U £ \ 2 dx. (1.2) 

z JT d xR d z J T d 
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In this paper we study the behavior of solutions to the system ED as e goes to 0. We will 
refer to this limit as the quasineutral limit. Let us notice that the Debye length is almost always 
very small compared to the typical observation length, so the quasineutral limit is relevant from 
the physical point of view and widely used in plasma physics. 

Let us observe that, if f £ —>• / and U £ —» U in some sense as e —» 0, the formal limit of our 
system is 

dtf + v • V*/ + E ■ V v f = 0, 


E = -V X U, 
fu d f dv = 1, 

k f\t=o = fo > o, / T d xR d fo dx dv = 1, 
and the total energy of the system reduces to the kinetic part of ED 


(1.3) 


£(/(*)) : =\[ f\v\ 2 dvdx. 

1 JT d xM d 


In this system, the force E is a Lagrange multiplier, or a pressure, associated to the constraint 
Jwj'fdv = I- 

The justification of the quasineutral limit from the rescaled Vlasov-Poisson system ED to 
ED is subtle and has a long history. Up to now, this limit is known to be true only in few cases 
and we refer to b 0 0 m m im e for a deeper understanding of this problem. 

One of the first mathematical works on the quasineutral limit of the Vlasov-Poisson system 
was performed by Grenier in [8]. He introduces an interpretation of the plasma as a superposition 
of a -possibly uncountable- collection of fluids and he shows that the quasineutral limit holds 
when the sequence of initial data /o j£ enjoys uniform analytic regularity with respect to the 
space variable. This convergence result has been improved by Brenier j3], who gives a rigorous 
justification of the quasineutral limit in the so called “cold electron” case, i.e. when the initial 
distribution /q j£ converges to a monokinetic profile 


fo(x,v') Po{%)d v = V Q(x) 

where 5 V denotes the Dirac measure in velocity. For further insight on this direction see also 

mmm- 

A different approach, more focused on the question of stability, or eventually instability, 
around homogeneous equilibria in the quasineutral limit is developed in E They show that 
the limit is true for homogeneous profiles that satisfy some monotonicity condition, together 
with a symmetry condition, i.e. when the initial distribution /o i£ converges to an homogeneous 
initial condition p(y) which is symmetric with respect to some v € M and which is first increasing 
then decreasing. 

In our previous work |12| . we consider the quasineutral limit of the one-dimensional Vlasov- 
Poisson equation for ions with massless thermalized electrons (considering that electrons move 
very fast and quasi-instantaneously reach their local thermodynamic equilibrium), and we prove 
that the limit holds for very small but rough perturbations of analytic data. In this context, small 
means small in the Wasserstein distance W \, which implies that highly oscillatory perturbations 
are for instance allowed. Our aim is to show that an analogue of this result holds in higher 
dimension. 


Let us introduce the notions of p -Wasserstein distance and weak convergence in the p- 
Wasserstein space V P {M) (see for instance (T9]h 
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Definition 1.1 (Wasserstein distance). Let (M,d) be a Polish space and let us denote with 
V P (M) the collection of all probability measures p on A4 with finite p moment: for some xq € M., 



d(x, xq) p dp(x) < +oo. 


Then the p-Wasserstein distance between two probability measures p and u in V p (A4) is defined 
as 


l/p 


W p (n,v):=[ inf d{x, y) p d^x, y) ) 

where T(p, v) denotes the collection of all measures on M. x M with marginals p and v on the 
first and second factors respectively. 


Definition 1.2 (Weak convergence in V p {M)). Let (A4,d) be a Polish space, and p € [l,oo). 
Let {pk)keN be a sequence of probability measures in V P (M) and let p be another measure in 
P(M). Then pk converges weakly in V p (M) to p if any one of the following equivalent properties 
is satisfied for some (and then any) xo € A4: 

1. ph —*■ p and f d(x,xo) p dpk(x) —>• f d(x, xfi) p d p{x)\ 

2. pk —^ p and lirnsup / d(x, xfi) p dpu{x) < j d(x, xo) p dp(x)\ 

k —^CXD 

3. p k p and lirn limsup f d(x )>R d(x, x 0 ) p dp k (x) = 0; 

it—>-°o k—too v — 

4. For all continuous functions ip with |<^( x ’)l — C(1 + d{x ,xq) p ), C € R, one has 

J y(x)dp k {x) ->• J ip(x)dp(x). 


In our case the Wasserstein space VziM.) is the space of probability measures which have a 
finite moment of order 2 and it will always be equipped with the quadratic Wasserstein distance 

W 2 . 

Remark 1.3. • By Holder’s inequality we have that 


p < q => W p < W q . 

In particular, the Wasserstein distance W\, is the weakest of all and results in W 2 distance 
are usually stronger than results in W\ distance. 

• Let (A4,d) be a Polish space, and p € [l,oo); then the Wasserstein distance W p metrizes 
the weak convergence in V p (Ai). In other words, if (pk)ke N a sequence of measures in 

Vp(M) and p is another measure in P(Ai), then pk converges weakly in V p (A4 ) to p if 
and only if 

Wp(pk,p)^ 0 as k^-oo. 

In order to state our main result, let us introduce the fluid point of view and the convergence 
result for uniformly analytic initial data introduced by Grenier in [8]. 

We first recall the definition of spaces of analytic functions we use below. 
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Definition 1.4. Given 5 > 0 and a function g : T —>• R, we define 


llslk ~ E l9(*)l« |t| . 

fcez 

where g(k) is the k-th Fourier coefficient of g. We define B$ as the space of functions g such 
that \\g\\B s < +oo. 

We assume that, for all e € (0,1), go j£ (x,v) is a continuous function; following Grenier [8], 
we write each initial condition as a “superposition of Dirac masses in velocity”: 



with M := R d , dp(6) = c d l+ d9 d+1 1 where c d is a normalizing constant (depending only on the 
dimension d), 

Po,e = -rt 1 + \°\ d+1 )goA x i v o,s = d • 

Of 

This leads to the study of the behavior as e —>• 0 for solutions to the multi-fluid pressureless 
Euler-Poisson system 

d t pi + V x ■ (p 9 E v 9 £ ) = 0, 
d t v 9 + v 9 ■ V x v 9 = E e , 

< E e = -V X U £1 (1.4) 

-s 2 A x U £ = f M p 9 dp(6)- 1, 

n 6 1 _ n 0 n .O | _ n .Q 

. Pe\t =0 Po,£i V £\t=Q v 0,s' 

One then checks that defining 

9e(t, x,v)= / p 9 (t, x)6 v=v eu } dp{6) 

JM 

provides a weak solution to (ED. 

The formal limit system, which is associated to the kinetic incompressible Euler system ED- 
is the following multi fluid incompressible Euler system: 

' d t p 9 + V* • (p 9 v 9 ) = 0, 
d t v 9 + v 9 ■ V x v 9 = E , 

< curl El = 0, f Td E dx = 0, (1.5) 

Jm P e d P( 0 ) = 

n 9\ _ _ oft 

P |t=o — Po,v k=o — ^0! 

where the p$ are defined as the limits of p^ £ (which are thus supposed to exist) and v 9 = 9. 

As before, one checks that defining 

g(t,x,v)= / p 9 (t,x)d v=v e {ttX) dp(0) 

J M 

gives a weak solution to the kinetic Euler incompressible system ED- 

We are now in position to state the results of [8, Theorems 1.1.2, 1.1.3 and Remark 1 p. 
369], 
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Proposition 1.5. Assume that there exist 5o,C,r] > 0, with p small enough, such that 


sup sup(l + v 2 )\\g 0t£ (-,v )\\ BSi < C, 
ee(o,i) ueR 


and that 


Denote for all 6 E M, 


sup 

e€(0,l) 


gO,e(-,v)d,V ~ 1 


< T], 


b 5 0 


Po, E = tt(1 + 6 2 )g 0t£ (x, 9), v 6 0e = v e = 6. 

Assume that for all 0£f, pg £ has a limit in the sense of distributions and denote 



Then there exist <5i > 0 and T > 0 such that: 


• for all e E (0,1), there is a unique solution (p £ , v £ )g e M of (11.41) with initial data (p(j e , v(j £ )geM> 
such that p £ ,v £ E C([0, T]; 13$^ for all 9 E M and e E (0,1), with bounds that are uniform 
in e; 

• there is a unique solution {p 6 ,v e )e^M of (jl.5[) with initial data (pg, Vq)q & m, such that 
p\ v e E C([0, T ]; B $ 1 ) for all 9 E M; 


• for all s E N, we have 


sup sup 
6£Mte[0,T] 



Q | | \\ 0 

P ll/f s (T) + | he 


1 it 

— (d+(i, x)e 

i 


d-(t, x)e ^)-v H \\ H s m 


0 


( 1 . 6 ) 

where d±(t,x ) are the correctors introduced to avoid the so called “plasma oscillations”. 
They are defined as the solution of 


curl d± = 0 , 


div d t d± + / p e v e p{d9 ) • V ) ri± ) = 0 


(1.7) 


where j £ 


div d±( 0 ) 

/ Pe v eP(d9). 


lim div 
£—^0 


^( 0 ) ± *J £ ( 0 ) 

2 


( 1 . 8 ) 


Remark 1.6. If in (ll. 8 [) . divd±(0) = 0, then the initial data are said to be well-prepared and 
there are no plasma oscillations in the limit e —>• 0 . 


The main result of this paper is the following: 

Theorem 1.7. Let 7 , <5o, and Cq be positive constants. Consider a sequence (fo, £ ) of non¬ 
negative initial data in L 1 for (ED such that for all £ E (0,1), and all x E T d , 

• (uniform estimates) 

||/o,e||oo < Co, £(fo,e)<Co, 
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(compact support in velocity) 


foA x ’ v ) = ° */ M > pf’ 

• (analytic + perturbation) There exists a function ip : (0,1] —> R + , with lim e _>.o t( £ ) = 0 
such that the following hold. Assume the following decomposition: 

fo,e = 90,e “I - ^0,e> 

where (go,e) is a sequence of continuous functions satisfying 

sup sup (1 + |u| 2 )||5o, £ (-,'y)||B5 n < C, 
ee(o,i) vm d 

admitting a limit go in the sense of distributions. Furthemore, (/io,e) is a sequence of 
functions satisfying for all e > 0 

W 2 (fo,e,90,e) = V( £ )- 


For all e € (0,1), consider f e (t) a global weak solution of (11,11) with initial condition /q. £ , in 
the sense of Arsenev m- Define the filtered distribution function 


fe(t,x,v) := f £ (t,x,v- -(d+(t,x)e^ - d-(t, x)e 


(1.9) 

\ i / 

where (d±) are defined in (11.71) . 

There exist T > 0 and g{t) a weak solution on [0, T] of (1 1.3 [) with initial condition g q such 

that 

lim sup Wi(f e (t),g(t)) = 0. 

£_> 0 te[o,T] 


Explicitly, we can take 

• in two dimensions, <p(e) = exp 


K 


eX P 1 £ 2(l+ma x(/3, 7 )) 


, for some constant K > 0, (3 > 2; 


• in three dimensions, (p(e) = exp 



K 

g.2+max(38,37)) 


for some constant K > 0. 


Remark 1.8. Let us notice that in Theorem m we consider sequences of initial conditions with 
compact support in velocity (yet, we allow the support to grow polynomially as s —>• 0). The 
reason is that, in the spirit of )12f . we rely on a Wasserstein stability estimate to control the 
difference between the unperturbed analytic solution and the perturbed one. In dimensions 2 and 
3, as we shall explain below, we need L°° bounds on the densities of both solutions. In order to 
have such a bound on the L°° norm of the densities we need to control the support in velocity. 
Such a condition was not required in our previous paper m since, in the ID case, we coidd use 
a “weak-strong” Wasserstein stability estimate and only a L°° bound on the imperturbed solution 
was needed. 


Remark 1.9. In the opposite direction, we recall that in the one dimensional case there is a 
negative result stating that an initial rate of convergence of the form ip(e) = £ s for any s > 0 is 
not sufficient to ensure the convergence for positive times. This is the consequence of instability 
mechanisms described in m and m, As a matter of fact, we expect that an analogue of this 
result holds also in higher dimension. 


6 








2 Overview of the paper 


The following of the paper is entirely devoted to the proof of Theorem 11.71 Let us describe the 
main steps that are needed to achieve this convergence result. 

1. We first revisit Loeper’s Wasserstein stability estimates m on the torus T d and with 
quasineutral scaling, which allows us to control W 2 (fi, fo), where /i and /2 are two given 
solutions of ed in terms of the initial distance W 2 (/i( 0 ), / 2 ( 0 )) and of the L°° norm of 
the densities p\ = f Rd fi dv and p2 = f K d f2 dv. This first step is performed in Section 13.11 

2 . In the one dimensional case studied in our previous work [T 2 ] we had a “weak-strong” 
type stability estimate; as a consequence a control of the L°° norm of the density of the 
perturbed solution f e (following the notations of Theorem II ,T[) was not required. 

In the higher dimensional case under study, such an estimate is needed. To achieve this, we 
give quantitative estimates of the growth of the support in velocity for solutions of ED- 
We separate the 2 and the 3-dimensional case since different tools are involved. 

While in the two dimensional case studied in Section 13.21 only elementary considerations 
are needed, in the three dimensional case, we shall use a more involved bootstrap argument 
due to Batt and Rein [2], see Section 1331 

3. We finally conclude in Section []] by combining the results of the two previous steps and 
Grenier’s convergence result stated in Proposition 11.51 


3 Proofs of Steps 1 and 2 


3.1 W 2 stability estimate 

We start by giving the relevant W 2 stability estimate, adapting from the work of Loeper [14] . 
Theorem 3.1. Let f \, f '2 be two weak solutions of the Vlasov-Poisson system ED, and set 


Define the function 



P2 



dv. 


Ait ) := 


l + “2 a/|I/ 3 2(0IIl°°(T'*) 


max {11 pi (t) 11 L oo (Td) , 11 p 2 (t) 11 L oo {T d)} 


(3.1) 


and assume that A(t) E L 1 ([0, T]) for some T > 0. Also, set 

F t [z\ : = 16 de log (i^W c °/o+ s H Vz E [0,d], i E [0,T], (3.2) 

Then there exists a dimensional constant Cq > 1 such that, if I+2(/i(0), / 2 (G)) < d, then for all 
t € [0, T], 


W 2 (/i(t),/ 2 (t))< 


F t [W 2 (fi(0), f 2 (0))] if F T [W2(fi(0),f 2 ( 0 ))] < d, 

d eP° Jo A( s ) ds ^ _p T [w 2 (/i(0), / 2 (0))] > d. 


(3.3) 
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Proof of Theorem 1 3.1\ Before starting the proof we recall two important estimates that follow 
immediately from [HI Theorem 2.7] and the analogue of |14t Lemma 3.1] on the torus (notice 
that \x — y\ < Vd for all x, y € T d ): 

Lemma 3.2. Let dq : T d —>• M solve 

-e?AVi = pi-1, * = 1 , 2 . 


Then 


1/2 


e 2 || Vv&h - Vd' 2 || L 2 (Td) < ^nax{||pi|| L oo (Td) , ||p 2 ||z,«>( T ^) }J W2(pi,p2), 

4 Vd 


£ 2 \Wi(x) - VT, : (y)| <C\x-y\ log ^_^ || Pi - l|| L °o (T d) Vx,y € T d , i = 1 , 2 . 
To prove Theorem 13.11 we define the quantity 

Q(t) := f 1 \ Yi(t, S\(s)) - Y 2 (t, S 2 (s ))| 2 ds 
Jo 

where 

Si, S 2 : [0,1 ] -> T d x R d 

are measurable maps such that ( Si)#ds = fi(0) and 

^ 2 (/i(0),/ 2 (0)) 2 = f 1 | Si(s) - S 2 (s)\ 2 ds, 

Jo 

while Yi = (Xj, Vi) solve the ODE 

Xi = V, 

Vi = -V*i(t,Xi) 

with the initial condition 1^(0, x,v) = (x,v). 

Thus, thanks to |14l Corollary 3.3] it follows that fi{t) = Yi(t)#fi( 0) = \Yi(t, Si)]#ds. 
Then we compute 

\jt Q{t) = Sl) ~ X2(t ’ 52)1 [Vl{t ’ 5l) " V5a)] dS 

+ f\vi(t,Si) -V 2 (t,S 2 )] [V*i(t,Xi(t,Si)) - V* 2 {t,X 2 (t,S 2 ))] ds. 

Jo 

By Cauchy-Schwarz inequality, we have 

Id 


2 dt 


Q(t) < \ / |X 1 (t,Si)-X 2 (t,S 2 )| 2 dsW / |l/i(L5 1 )-E 2 (L5 2 )| 2 ds 


+ 


^ |Vi(i, Si) - V 2 (t, S 2 )| 2 ds^£ |VTi(t, X x (t, Si)) - X* 2 (t, X 2 (t, S 2 ))\ 2 ds 


< Q(t) + VoU)\ / -V^i(i,X 2 (t,5 2 ))| 2 ds 


| VTi (t, X 2 (t, S 2 )) - V ^ 2 (t, X 2 (t, S 2 )) | 2 ds. 


+ Vo(t) 













Using the definition of the push-forward, we finally get 


1 d 

2 dt 


Q(t) < Q(t) + \/Q(*)yj[ -VV 1 (t,X 2 (t,S 2 ))fds 


/T d xR d 


|V'Ui(i,rc) — x)| 2 / 2 (t, x,u) dxdv 


< Q(t) + VQ^)^jJ o -v^ l (t,x 2 {t,s 2 ))\ 2 ds 

+ ^\\p 2 (t)\\ L 0 0{T d ) ^/Q{t)^l |V^i(t,x) - X^> 2 (t,x)\ 2 dx. 


We now apply Lemma [3.21 to the last two terms and we bound them respectively by 

4 yfd 


C 


and 


\\Pi(t) lh-w jf 1 | Xl (f, 5l) _ x 2{t , 5 2 )| 2 log 2 


S 1 )-X 2 (t,S2)\ 


ds 


IM^Hiz-crd) max{||pi(t)|| L oo (T d), H^WIIl-^)} W 2 (j>i(t), p 2 (t)) 


1/2 


Since W 2 {pi(t), p 2 (t)) < \/Q(i) (see for instance ini Lemma 3.6]) we conclude that 




+ C 


1 + ^2y IIP2(0llL°°(T d ) m ax{ ||pi (t) ||l / 00 (T ti ) i I|P2 (i) ||L°°(T d ) } 

H^r(^) — ^.2 II J ' 0 ° (T d ) ^i) ~ ^ 2 >1 2 log 2 ( 


1 / 2 ' 


Q{t) 


4 y/d 


,S 1 )-X 2 (t,S 2 )\ 

Noticing that |Xi(t, Si) — X 2 (t, S 2 )\ < Vd (since X\ and X 2 are points on the torus) and 


ds. 


, /4 Vd\ 1 /16d 

lo g( — = o lo S 


V z > 0, 


we get 


£\X 1 (t,S 1 )-X 2 (t,S 2 ) | 2 log 2 (^- 


4\/d 


ds 


\Xi(t, Si) — X 2 (t, 5 2 )|“ log 


5i)-X 2 (/,5 2 )|. 

2 , 2 / 16 d 


|W 1 (L5i)-X 2 (t,5 2 )p 
1 


ds 


4 y„ 9(s) log 1 *' 


16 d 


where we set g(s) := |Xi(t, Si) — X 2 (t, -S l 2 )|" 
Hence, since the function 


= J " 1 °g 2 ( 1 w) ^0<z<d, 

| d log 2 (16) for z > d, 


(3.4) 
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is concave and increasing, recalling that g < d and applying Jensen’s inequality to H we get 


* 


1 


1 + ^2 y ||P2(i)||l,oo(T d ) max{ ||pi(f)|| L oo( T d), \\P2(t) || L oo( T d) } 

+ C 


1 / 2 ' 


Q(t) 


l,Pl(t) £ l h °° {Td) VQ{i)^H^f\( S )ds^ 


< 1 + 


1 1/2 


IIP2 (*)||l<»(T‘') m ax{||pi(t)|| L oo (Tti) , 11/92(i)llL-(Td)} )Q(t) 


+ 




where for the last inequality we used that g(s ) ds < Q(t). 

In particular, recalling the definition of A(t) in (13.ID . by (13.411 we deduce that there exists a 
dimensional constant Cq > 0 such that 


■^Q(t) < C 0 A(t)Q(t) log 


16 d \ 

I as long as Q(t) < d, 


while 


—Q(t) < Co A(t) Q(t ) when Q(t) > d. 


I 11 particular, assuming Q(0) < d, by (13.51) we get 

Q{t) < l6de 1 ° g (^) exp [ C ' 0 ^ A(s) ^ =: F t [Q(0)] 


(3.5) 

(3.6) 

(3.7) 


as long as Q(t) < d, which is the case in particular if Ft[Q{ 0)] < d. On the other hand, if there 
is some time to such that F to [Q( 0)] = d, since Q(to) < F to [<3(0)] by (13.61) we get 


Q(t) < d e Co J*o A(s) ds < d e Co A(s) ds for t > t 0 . 


(3.8) 


Noticing that Ft is monotone in t, we deduce in particular that if i^[Q(0)] < d and Q(0) < d 
then (13.71) holds, while if i^[Q(0)] > d and Q(0) < d then one can simply apply (13.81) . Finally, 
if Q(0) > d then we apply (13.61) to get 

Q(t) < Q(0)e c °ti A ( s '> ds . 

Combining these three estimates and recalling that Q(0) = IF2 (/i(0 ), / 2 (G )) 2 while Q(t) > 
this concludes the proof. □ 


3.2 Control of the growth of the support in velocity in 2D 

In this section, d = 2. Our goal is to obtain estimates on the growth in time of the support in 
velocity. This allows us to get bounds for the L°° norms of the local densities on some interval of 
time [0, T\. Recall that in the end, they will be used to apply the Wasserstein stability estimates 
proved in Section 13.11 

For f £ a solution of m, define 

V e (t) := sup {|n| : v € M 2 , 3x € T 2 , f e (t , x, v) > 0} . 

The key point of this Section is the following Proposition. 
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Proposition 3.3. Suppose that 


||/e(0)||oo < Co, J (\v\ 2 + U e (0, x)^f £ (0, x,v)dvdx<C 0 . 

Assume that 14(0) < Co/e 1 , for some 7 > 0. Let T > 0 be fixed. For all ft > 2, there is Cp > 0, 
such that we have for all e E (0,1) and all t E [0, T\, 

V £ {t) < (3.9) 

Therefore, for all fi > 2, there is C^ > 0, such that, for all e E (0,1) and all t E [0, T], 

llPelloo < C'p/e 2max Wrt. (3.10) 

In order to prove this Proposition, we shall use for convenience the change of variables 
(t,x,v) i->- (-, - ,v ). This leads us to consider, the following Vlasov-Poisson system, for {x,v) E 


±T 2 x M 2 : 
£ 


dt9e T V ' V xSe T l ; • V v9e — 0, 

F e = -V x $s, 

= Jig2 9e dv flj 2 x jj2 ffe dv dx , 

£ 

g £ \t .=0 = 90,e > 0, 

,/it2 x r 2 9o,e dx dv = -pz . 


(3.11) 


We shall denote 


and define for all t > 0 , 


■-S 


Ve ■= / 9e dv 


1 . 


V{t) := sup < |u|, v E Kr, 3x E -T 2 , g £ {t, x, v) > 0 > =V £ [ - 


(3.12) 


In the following, for brevity, the notation p , for p E [l,+oo], will stand for L p (^T 2 x ]R 2 ) or 
L p ( 7 T 2 ), depending on the context. 

The main goal is now to prove the following Proposition, from which it is straightforward to 
deduce Proposition 13.31 by applying the result for t = —. 


Proposition 3.4. Let /3 > 0. Let Cq > 0 such that for all e E (0,1), 


||Se(0)||°° < Co, f (\v\ 2 + $ £ (0,xj)g e (0)dvdx 

I It 2 xi 2 ' 7 


Co 

~ c -2 ’ 


< 


(3.13) 


Then for all a E (0,1), there exist C a , C' a > 0 such that for all e E (0,1), and all t > 0, 


V(t)< 


c n 


rQ +1 


t+ii+vm 


1 —a 


i/(i— a ) 


— 1 . 


(3.14) 


We will use the following standard property of conservation of L p norms and energy for 
solutions to the Vlasov-Poisson system in the sense of Arsenev: 


Lemma 3.5. For all t > 0, we have 
l|fi , e((')l | 00 < Co, I 


-T 2 xR 2 


^|u| 2 -I- $ e (t,x)jg e (i) dvdx < 


(3.15) 
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We also rely on the following Lemma about the Green kernel of the Laplacian on ^T 2 
(obtained from standard results on the Green kernel of the Laplacian on T 2 , after rescaling). 
We refer for instance to Caglioti and Marchioro [5ij. 

Lemma 3.6. There exists I\$ € C'°°(T 2 ;M 2 ) such that, denoting 

1 X 

K e( x ) = +£ K o(ex), 

we have for all x € [—1/e, 1/e] 2 , 

F e {x) = [ K e (x - y) \j] e {y) ~ l] dy. 

The key ingredient is the following Lemma, in which we obtain some appropriate 00 bound 
for the electric field, which allows us to control the growth of the support in velocity. 

Lemma 3.7. We have the following bounds. 

1. There is a constant C\ > 0 such that for all e € (0,1), and all t > 0, 

he{t)h<^, hem-KC^it) 2 . (3.16) 

2. There is a constant C 2 > 0 such that for all e € (0,1), and all t > 0, 

1/2 \ 

j • (3-17) 

3. For any a > 0, there is a constant C a such that for all 0 < t' < t, 

V(t) < V{t') + f (1 + V{s)T ds. (3.18) 

£ Jt> 

Proof of Lemma \ 3 . 71 In this proof, C > 0 will stand for an universal constant that may change 
from line to line. 


I FS,‘ 


< C 2 1 + - 


log-(l + F(i)) 


1. By the following interpolation argument, we have for all R > 0 


Ve = 




g £ dv + 



9edv < || g £ 


,R 2 



v\ 2 g £ dv, 


so by optimizing with respect to R we deduce that there is a C > 0 such that for all 

1 tt2 


t > 0,3:6 ^T 2 , 


\Ve\(t,x) < C g £ (t,x,v)\v\ 2 dv^j 
By (13.1511 . we deduce the first estimate of (13.161) . 


\ 1/2 


For what concerns the 00 estimate for g £ , it is a plain consequence of the definition of V(t), 
which controls the support in velocity. 
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2. By Lemma ETUI there holds for all t > 0, x € [—1/e, 1/e] 2 


I F s \(t,x)<— — 

2 vr J[- 1/£y i / E ]2 |x - a:' 


-|r? £ (t,a:') - 1| dx' + e||/\ 0 ||oo||(% - l)||i 


We observe that, since ||(% - l)||i ( i T2 ) < ||%||i(i T 2 ) + ||l||i(i T 2 ) = 

£||A'o|U||fe-l)||. <§. 

Let R > 0 to be fixed later. We have, using the Cauchy-Schwarz inequality, 

r i 


'[— 1/e,1/e] 2 I"*' ® 


7 || Ve(t,x') - l|dx' 

^— r \rj E (t,x')-l\dx'+ [ - — ^-fi\ri e (t,x')-l\dx' 

-* J\x-x'\>R\ x - x \ 


J\x—x’\<R 1 ^ 

< CR(\\r) e \\°° + 1) + ( 11%|| 




< 


Ci?(||%||o° + 1) + ( 





1/2 


„/|2 


dx' 


We choose R = „ ,/ +1 , which yields, 


1 



\\F £ \\oo < C 1 + ||%||2 + - log -(1 + ||??||°°) 


Using Point 1., we obtain the claimed estimate. 

3. Let a > 0. Let 0 < t' < t and let (x,v) such that g £ (t',x,v) ^ 0. Introduce the 
characteristics (X(s, t 1 , x, v), £(s, s, x, v)) satisfying for s > t' for the system of ODEs 


—X(s,t',x,v) = £(s, t', x, v), X(t',t',x,v) = x, 

as 

—£(s,t',x,v) = F £ (s,X(s,t',x,v)), £(t',t',x,v) = v. 

as 


(3.19) 


We have 


£(t, t', x, v) = v + f F £ (s,X(s,t',x,v)) ds. 
Jt' 

Therefore, we have, using Point 2., 

If I (t,t',x,v) < |u| + J \F e \(s,X(s,t',x,v))ds 


< \v\ + — f | 1 + 

£ Jt’ 


log-(l + U(s)) 


1 / 2 " 


ds. 


Since g £ satisfies (13.111) . it is thus constant along the characteristics (13.191) . and we have 

g £ (t, X(t, t', x, v),£{t, t',x, v)) = g e (t',x,v). 
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By definition of V[t) and V(t'), we obtain 


v(t) < v(t') + y /* (i + 


log-(l + V(s)) 


1 / 2 " 




In order to get the polynomial bound on V(t), let r a > 0 such that for all x > 1, 

(log x ) 1,/2 < 1 + r a x a . 


We thus get 

V(t)< V(t') + q ’ (2 1 + r ° ) /‘(l + V(»))“do, 
e Jt' 

which proves our claim, taking C a := 6 * 2(2 + r Q ). 


□ 


Equipped with this result, we can finally proceed with the proof of Proposition 13.41 

Proof of Proposition \3.J\ We begin by observing that dividing by l/(t — t') in both sides of 
(13.1811 and letting t' —>• t we deduce that 

fv(t) £ + r(«))“. (3.20) 

We will obtain the claimed bound by a comparison principle. To this end, introduce a small 
parameter p, > 0 and define 

wwt) ~ ( 2[C ° ~ a) t + (1 + V(0) + m) 1 -”) ■ 

By construction, it satisfies for t > 0 

|«W«) = (3.2i) 

and 

W M (0) = 1 + E(0) + fr > 1 + 1/(0). 

We claim that W fl (t) > 1 + V(t) for all t. Indeed, let 

t 0 := inf{t > 0 : W,,_({) < 1 + V(t)}, 

and assume by contradiction that to < +oo. Notice that because /i > 0 we have to > 0. Then 
by continuity at the time to we get 


W M (t 0 ) = 1 + V(t 0 ). 


By (13.201) and (13.211) . we have 

d 


dt 


(V(t) - Wfj,(t) - l) lt=t0 < -g^(l + V(t 0 )) a - %^^(t 0 ) c 




— pl+a 


W M (t 0 ) Q < 0. 
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This is a contradiction with the definition of to- 

Hence we obtained that for all fi > 0 and all t > 0, 

1 + V(t)<W /i (t). 

By taking the limit p, —>• 0, one hnally gets for all all t > 0 

1 + v «) £ ( 2C °l+7 + (i + rfo)) 1 -”), 

which proves the Proposition. □ 


3.3 Control of the growth of the support in velocity in 3D using Batt and 
Rein’s estimates 

In this section, we deal with the case d = 3. We consider as before, for f e a solution of (D>- 
V £ {t) := sup{|n|, v G K 3 , 3x G T 3 , f e (t, x, v) > 0} . 

We have in 3D the analogue of the key Proposition 13.31 in 2D. 


Proposition 3.8. Suppose that 

11/e (0)1 loo < Co, J (|u| 2 + U E (0,xjjf E (0,x,v)dvdx < Cq. 


Assume that 14(0) < Cq/e 1 , for some 7 > 0. Let T > 0 be fixed. There is C' > 0, such that we 
have for all £ G (0,1) and all t G [0,T], 


V e {t) < 


C' 

gmax{38/3,7} ' 


(3.22) 


Then 


11 Pe 11 OO 4 


C' 

g-max{38,37} ’ 


(3.23) 


Note that this result involves exponents which are “more degenerate” than in the 2-D case. 
We shall prove this result as an application of the estimates obtained by Batt and Rein in |2j. 
The result of |2j is an adaptation to the case of the torus T 3 of the fundamental contribution 
of Pfaffelmoser m (see also [181 T3]h which allowed to build global classical solutions of the 
Vlasov-Poisson system in M 3 xl 3 . In M 3 xl 3 , it may be possible to get better estimates than 
(13.221) (i.e. with smaller exponents) by using dispersive effects, see [171 H~5l 13j and more recently 

m- 

In order to prove this Proposition, we shall use the change of variables (t,x,v) > (|, |,u). 

This leads us to consider, the following Vlasov-Poisson system, for (x,v) G ^T 3 x M 3 : 


dt9e T V ■ V x9e T F e ■ S7v9e — 0) 

F e = -V x $ e , 

= Jjj3 9e dv — fij3 x jj3 9e dv dx , 

£ 

9e\t=o = 9o,e F 0, /i^3 x r3 9o,e dx dv = -^ 5 . 


(3.24) 
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We shall denote as in the 2D case 


Ve ■= / 9e dv, 
J R 3 


and define for all t > 0 , 


V(t) := sup 


|u|, v £ M 3 , 3x € -T 3 , g e (t, x, v) > 0 


(3.25) 


As before, in what follows we use the notation p , for p £ [1, +oo], will stand for L p (^T 3 x M 3 ) 
or L p (^T 3 ), depending on the context. 

The main goal is now to prove the following Proposition, from which we deduce Proposition 
13.81 by choosing t = j. 


Proposition 3.9. Suppose that 


||/ £ (0)|h <C 0 , 



e 


Let 7 > 0. Let Cq > 0 such that for all e £ (0,1), 


1 /( 0 ) < 


Co 

e'y ' 


Then there exists C\ > 0 that for all e € (0,1), and all t £ [0, T\, 

Ci 


V(t) < rnax{^ + 


r32/3 


+ ,LL t . + 4 - 5 l _ 

V e 64 / 3 g-32/3+7 


^0 +r -7/2l 

e'y V 


(3.26) 


(3.27) 


Proof of Proposition \S.9\ Consider the usual notations for characteristics of (13.1911 and intro¬ 
duce, as in [ 2 ], 


hi(t) ■.= sup{||%(s)||loc, 0 < s < t} + 1, 

/ 12 (f) := sup{|^(s, t,x,v) — v |, 0 < s,r < t, (x , v ) £ -T 3 X M 3 }. 

We look for a bound on / 12 , which will imply the control on V(t). To this end, we crucially rely 
on the key bootstrap result in the paper of Batt and Rein |2j, which we recall in the form of a 
lemma for the reader’s convenience. 

Lemma 3.10 (Batt, Rein). Assume that there is C* > 0 and /3 > 0 such that 

/ 12 (f) < C*thi(t)P 

then for some universal constant C > 0 (that hereafter may change from line to line), 

h 2 {t) < Cf(c* 4 / 3 /if /3 (f) + 1 (V 6 (t) + ), (3.28) 

if hi(ty 13/2 < t. 
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By using [21 Eq. (5), Section 4, p.414], there is C > 0 independent of e such that for all 
e > 0 and t > 0, 

h 2 (t) < Cthi(t) 4/9 . (3.29) 

We deduce from Lemma [3.101 and fl3.28|> that 

h 2 (t) < ~^thi(t) 8 / 27 , if hi(t) -2 / 9 < t. 

E 6 

Using fl3.28|> twice, we finally obtain 

h 2 (t) < ^thi(t) 16 / 81 , if /ii(f) -4 / 27 < t 

and since ^ < |, we get 

h 2 (t) < if hi (t)~ 8/81 < t, 

Using the straightforward bound 


hi{t) < C (U(0) + h 2 (t)) 3 , 
we deduce that for all e E (0,1) and t > 0, 




1/2 


if h\{t) 8 / 81 < t. On the other hand, if hi(t) 8 / 81 > t then, 


M*) < ^sTsT 


and thus, by (13.2911 . we get 


h 2 {t) < Ct~ 7/2 . 

We conclude that for all e E (0,1) and t > 0, 


f 1 / 1 \ 1/2 . 

h 2 (t) < Cmaxj-^f + h 2 {t)J , t' 7/2 |, 


which yields, taking t = T, 


h 2 (T ) < max / - 


Therefore, for iE [0, T] we get 


V(t) < max | ^ + 


C ’ + J^Lt‘ + 4- °' 


g-32/3 ' y £ 64/3 1 g-32/3+7 


T~ 7 / 2 L . 


C ‘ + . TLt* + 4 c ‘ 


r 32/3 ' 1/ p64/3" 


£ 32 / 3+7 


O) 

’ e 7 


+ T - 7 / 2 , 


which proves Proposition 13.9 


□ 


17 


























4 Proof of Theorem 11.71 


We prove the main Theorem by a perturbation argument, relying on the Wasserstein stability 
estimates of Theorem 13.11 

Before that, we first state a Lemma about the effect of ^-dependent translations in the 
velocity variable on the W\ distance. 

Lemma 4.1. Let /j, v be probability measures on Vi(T d x M . d ) and let ft, D be probability measures 
on V\(T d X W l ) defined as follows: 

< ft, (j)(x, v) >:=< ft, cj)(x, v — C(x )) > for all <f> € Lip (T d x M rf ); 

< v, fi(x, v) >:=< v, cf(x, v — C(x)) > for all cf> € Lip (T d x M d ). 

Then 

Wi(ft, u) < (1 + \\D x C\\ l °°) Wi(ft,v) where D X C := {p Xi Cj{x)) 0<i j<d . (4.1) 
Proof. By the Kantorovich duality we have the following expression: 

Wi(ft,D)= sup 

Let us denote 0(£,u) = (p(x,v + C(x)) and computing the gradient of 0 we deduce that 

HV’lliip — (l + \\D x C\\l°o) \\tp\\Lip < (l + \\D x C\\ L oo) 

from which we obtain (14.11) . □ 

We can now proceed with the proof of Theorem 11.71 Let /o,e,5o,ei ho,e satisfy the hypotheses 
of Theorem 11.71 Using the same notations of the statement, we want to show that for some 
T > 0, 

lim sup Wi(f e (t),g(t)) = 0. 

£_! ’° te[o,T] 

In analogy with the definition of f e we define 

g e (t,x,v)= / p e £ (t,x)5 v=v eu x) dg{6) 

Jm 

and 

9 s(t, x,v)= / p d £ {t, x)5 v=v e {ttX)+Ce ^ x) dp(0) 

J M 

i it _ it 

where C £ {t,x) := — i(d+(t, x)e V® — d-(t,x)e V®). 

We now prove the following estimate: 

WfiJ^g) < Wi(f £ , g e ) + Wi(g £ ,g). (4.2) 

We first consider the second term in the right hand side. The uniform convergence to 0 follows 
from Proposition 11.51 and the Sobolev embedding theorem. We have indeed for some T > 0, for 


J g>{x, v + C(x))dp(x) — g>(x, v + C(x))d v(x) 


18 






all t E [0, T]: 


Wi{g e (t),g(t)) = sup (g £ - g , <p) 
lbl|Lip<l 

= sup l / (pf(i, x)p(x, v £ (t, x) + C £ (t, x)) — p e (t, x)<p(x, v e (t, x))) dp(0) dx 

IHlLi P <i Ut*Jm 

= sup < / / p £ (t, x)(<p(x, v £ (t, x) + C e (t, x)) — p(x, v e (t, x)) dp{0) dx > 

llvllLip<l UT d JM " J 

+ sup l / (p £ (t,x) - p 9 (t, x))ip(x,v e (t,x)) dp(0) dx\ . 

IM|u P <l wT d JM J 


Thus, we deduce the estimate 

Wi(g e (t),g(t)) < sup sup \\p e £ \\oo\\v\\u P [ \\v £ (t,x) + C £ (t,x) - v e (t,x)\\ oo dp(0) 
IM| L i P <i e 6 (o,i), e&M Jm 

+ sup / \\p £ - pe\\oo dp(0)\\p\\u P ( 1/2 + sup ||/(t,x)|| 00 ) 

NlLi P <i^A4 V e&M J 

+ sup | j j (p £ {t, x) — p e (t, x))<p(O,O)dp(0)dx 1 . 

IMlLi P <i Ur d xRJM " J 


We notice that the last term is equal to 0 since for all t > 0, 




p £ (t, x) dp{0) dx 



dx = 1, 


by conservation of the total mass. Considering the supremum in time, we see that the other two 
terms converge to 0 , using (USD- so that we get 


lim sup Wi(g e ,g) = 0. 
e_s ' 0 te[o,T] 

We thus focus on the first term of the right hand side of (14.21) . First we use Lemma 14.11 to see 
that 

Wi(/ e ,&) < (1 + \\D x C e (t,x)\\ L <x) W\(f £ ,g £ ). 

Observe from the definition of the corrector C £ , there is Ct > 0 independent from e such that 
for all t E [0, T], 

II D x C e (t, -)||z,°° < Ct- 

We therefore have to study W\ (f £ , g £ ). We first use the rough bound 


W 1 (f £ ,g £ )<W 2 (f £ ,g £ ), 
then use Theorem 13.11 to get the estimate 


sup W 2 (f £ ,g £ ) <16 d exp {log 
te[o,T] 

where 


16 d J 


exp 


1 


Co T ~2 ( 1 + H/ 5 /elU oo ([0,T];L“>) + IIPsJl°°([ 0,T];.L° o ) 


Pfe ■= / fedv, p 9e = / p e £ dp(0). 
jR d Jm 

Recalling (USD- we have for some C > 0 independent of e that 


\\Pg s \\L°°({o,T]-,L™) < C. 


For what concerns pf e we apply 
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in two dimensions, (13.101) in Proposition 13.31 to infer that for all j3 > 2, there is some 
C/s > 0 independent of e such that 


II Ps 


< 


Cp 


£ 2 max{/J, 7 > ’ 


in three dimensions, (I3.23[) in Proposition 13.81 to infer that for some C > 0 independent of 
e such that 


II Pe 


< 


C 


gmax{38,37} ' 


We deduce that choosing 

• in two dimensions, 99 (e) = exp 

• in three dimensions, 99 (e) = exp 


K 


exp --5(7+ 

max(/3,7)) 


eX P l e 2+max(38,37)) 


, for some constant K > 0; 

, for some constant K > 0, 


up to take a smaller time interval of convergence [0, T], 


We conclude that 


lim sup W 2 (fe,g e ) = 0 . 
£_> ' 0 te[o,T] 


lim sup Wi(f £ ,g) = 0 
£_) ' CI te[o,r] 


and the proof of Theorem 11.71 is complete. 
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